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1. Motivation 

Experiments with classical waves have become a very versatile tool to study localization 
due to disorder. In particular the experiments by Lagendijk and coworkers (Wiersma 
et al 1997) on the localization of light in powders, and by the Genack group on the 
localization of microwaves in disordered metallic spheres (Chabanov and Genack 2001) 
have to be mentioned (for a review of these types of experiments see Soukouhs 1996). 
Moreover microwave techniques are able to study spatially resolved field distributions 
in disordered systems of linear dimensions in the order of some 10 cm (KudroUi et al 
1995, Stockmann et al 2001). Such quantities are inaccessible in electronic quantum 
dot systems of submicron size (except for the recent experiments by Topinka et al 
(2001)). With increasing frequency one typically observes a transition from localized 
to delocalized wave functions, depending on the number of scatterers and the strength 
of the scattering potential. Pulse propagation can be studied as well by microwave 
techniques as has been shown by Stein et al (1995). All quantities of interest are thus 
experimentally accessible in disordered systems, including conductivity, localization- 
delocalization transitions, pulse propagation, transition from the ballistic to the diffusive 
regime, and so on. 

On the theoretical side the situation is less favourable. Though there is a vast 
amount of literature on disordered systems already in the seventies and the eighties of 
the last century (see e.g. Anderson 1978, Lee and Ramakrishnan 1985 for reviews), 
there is as yet no theory available covering the complete range from the localized to 
the delocahzed regime. Today the standard approach to study disordered systems uses 
supersymmetry techniques to arrive at Efetov's non-hnear a model (Efetov 1983). It 
has the serious draw-back that the occurring supersymmetric variables are field variables 
depending on the position. Only in the zero mode approximation, where the position 
dependencies are neglected, the model can be solved exactly and reproduces random 
matrix theory. This is why localization-delocalization transitions cannot be obtained in 
this way. Only perturbational corrections are possible, with the consequence that e. g. 
the distribution of wave function intensities deviates slightly from the Porter-Thomas 
behaviour found in the delocalized regime (see Guhr et al 1998, Mirlin 2000 for reviews). 

In this paper an alternative approach is proposed which avoids the complication of 
position-dependent supersymmetry fields. Moreover, it is even closer to the situation 
met in experiments, as Efetov's ansatz. 

2. The model 

Let us consider a biUiard system with hard walls and statistically distributed scatterers 
described by the Hamiltonian 



H^Ho + V, 



(2.1) 
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where Hq is the operator of kinetic energy with matrix elements 

(Ho)^^ = EXm (2.2) 

and V is the potential energy of the scatterers. In Efetov's approach the potential is 
assumed to be delta correlated, 

{V{r)V{r')) 6{r -r') , (2.3) 

which gives rise to the above mentioned problems with position-dependent 
supersymmetry fields. In this paper the more explicit ansatz 

L 

V{r) = ATiXY^Sir -n) (2.4) 

1=0 

is used instead, where the ri are the positions of the scatterers, and L is its number. The 
factor An has been introduced for later convenience. This ansatz dates back to Lifshitz 
(1964) and has been applied since then by various authors, among others Luttinger et 
al (1983, 1987). 

In the basis of eigenfunctions ipn{f) of the billiard without scatterers the matrix 
elements of V{r) read 

L 

K,n = 47rA5]C(n)^m(n). (2.5) 

To simplify the calculations we assume that time-reversal symmetry is broken, e. g. 
due to the presence of an applied magnetic field. 

We are now going to calculate the density of states, averaged over the positions of 
the impurities, 



(p(^)) = -^Im^Tr(^-^ 



(2.6) 



as 



HJ / ' 

where = E + ie. Using standard transformations, equation (|2.6[) can be written 



PiE) = -l^lm[Z{E,E')] 



where 



Z{E,E') 





-H 


\E-^ 


-H\ 



(2.7) 

E'=E 



:2.8) 



Z may be written in terms of an integral over commuting and anti-commuting 
variables as (Verbaarschot et al 1985) 

Z=(^Jd[x] exp (^i Y: [{E+6^^ - H^p) xlxp + {E'^5^p - H^p) C^/?] j ^ , (2.9) 
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where 



N 



= n dx*dxQd^*d^Q 



(2.10) 



a=l 



We adopt the usual convention and use Latin letters for commutating and Greek 



ones for anti-commutating variables. In short-hand matrix notation equation (|2.9| ) may 
be written as 



ixt(E-H)x 



where 



and 



d[x] e 



/ E ■■■ \ 

E = e®1n= ; ••. ; 

V ■■■ E ) 

( Hnl ■■■ H^nI \ 



E 



E, 



EL 



(2.11) 

(2.12) 
(2.13) 



H= 1®H 



(2.14) 



V H]yil ■ ■ ■ Hnn'^ J 



In Equations (p.l3|) and ( p.l4| ) l^v and 1 denote the A^- and the two-dimensional 
unit matrix, respectively. Inserting expression ( |2.1| ) for H, equation (|2.11|) reads 



d[x] e 



ixt(E-Ho)x 



-47riA J2 '4>^in)i'l3in){x''^xp+(^^0) ^ 



lal3 



where 



J d[x] e-^(E-Ho)x^L ^ 



/ -47TiXj2rjr)ijfi{r){x*^xp+^l^^)^ 

M = (e 



(2.15) 



(2.16) 



The average in equation ( |2.16|) has to be taken over the positions of the impurities. 
But, and this is the central idea of this paper, instead of varying over the positions, we 
may equally well calculate this average by weighting the expression on the right hand 
side of the equation with the joint probability density p^ipiRy'ipii, ■ ■ ■) to find at any 
point in the billiard the values ipiR, ipii^ • • • for the wave function amplitudes. If the 
billiard without scatterers is chaotic, the probability density factorizes, p{ipiR, ipn, ■ ■ ■) = 
Yipii^aRjpii^ai), and real and imaginary part of the wave functions are Gaussian 

a 

distributed, 



p{ipi 



lA 
— e 

71 



p(^/) 



TC 



(2.17) 
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where A is the bilhard area. The average ( p.l6| ) over the impurity positions may hence 
be written as 



M 



-47riA^ 1/.* V/3(a;* 

YlldiJaRdiJalPii^aRjpii^al)] e 



(2.18) 



With the weight functions ( |2.17|) the integrations are easily performed yielding 

1 



M 



where X is the N x N matrix with the elements 



(2.19) 



(2.20) 



According to the Weyl formula the mean density of states in two-dimensional 
billiards is given by (p) = A/Att. Following the usual practice we normalize this quantity 
to one, and omit the factor Atc/A in the following. The determinant (|2.19|) is now 
transformed by means of the relation 



\1n + AB\ 



BA\ 



(2.21) 



holding for arbitrary N x M matrices A, and M x N matrices B. This follows in a 
straightforward manner from the relation \M\ = exp[Tr(lnM)]. It is not necessary that 
the matrices are quadratic providing us with an efficient tool to reduce the rank of 
determinants. Applied to equation ( p.l9|) relation (|2.21| ) yields 

1 



M 



iXX 



(2.22) 



where X is the 2x2 supermatrix 



X 



We have thus arrived at the intermediate result 



d[x] e 



ixt(E-Ho)x 



1 + iAX 



(2.23) 



(2.24) 



Whenever there are super matrices involved, determinants and traces have to be 
interpreted as super determinants and super traces, respectively, where we shall use the 
convention of Verbaarschot et al 1985. 



It is instructive to consider the small A limit of expression (|2.24|) . The determinant 
may be expanded as 



1 + iAX 



exp 
exp 



-LTrln(l + iAX 
.iLATrX-^ 



TrX^ + . 



(2.25) 
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Stopping at the quadratic term, equation ( p.24| ) reads 



d[x] e 



ixt(E-Ho-LAljv)x--^ TrX 



(2.26) 



where TrX = x^x was used. This is exactly the expression obtained for the ensemble 
average of the Hamiltonian 



H = Ho + LAljv + Hi 



(2.27) 



where the matrix elements of Hi are Gaussian random variables with variance (Hf) = 
. We thus can note already at this early stage that in the small A limit random 
matrix results will be recovered. 



3. The X integrations 



The usual approach to perform integrations of type (|2.26| ) is a Hubbard-Stratonovich 
transformation to remove the TrX^ term in the exponent, depending on the integration 
variables in the fourth order. As a result the x integrations reduce to simple Fresnel 
integrals which are trivially solved. 

For the integral ( |2.24D a Hubbard-Stratonovich transformation is not possible. One 
way to proceed further is to write again the determinant in terms of a super integral, 

1 /■ ,r , -yt(i+iAx)y 



1 + iAX 



(3.1) 



We need L replicas of this integral since the determinant enters in the Lth power, 
leading to the introduction of 4L new integration variables. The x integrations can then 
be performed in the usual manner. 

To avoid the introduction of such a large number of new integration variables, we 
apply another approach. Let us consider the integral 



/z.(A)=|d[t]|T|VT^(^^), 



where 



A 



a a 
a a 



T 



t T* 

T i 



(3.2) 



(3.3) 



are supermatrices of rank 2. 

Because of the basis independence of trace and determinant it is always possible to 
chose the T variables such that A is diagonal, i. e. a = a* = 0. Equation ( |3.2| ) then 
reads 

L 



h{A) 



dt dt dr* dr 



t T* 
T t 



(3.4) 



Introducing new integration variables s = at, s = at, a* = ar* , a = ar we obtain 
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h{A) 



where 



- — - — \a\ da \a\ da 

\a\ \a\ 

\-\\ ^ 

n) 

a / 



d[t]|T|^e-'r^(^) 



a ^s 



a a 



-1 * 
a a 

d^^s 



- Tr{s-s) 



(3.5) 



or 



h{A)/h 



(3.6) 



(3.7) 



Because of the basis independence of this expression the latter result holds for 
arbitrary supermatrices A, not necessarily diagonal. 

This is an alternative to equation (|3.1| ) to express the power of a determinant as 
a superintegral, avoiding the need to introduce L replicas. The question arises, how 
the paths of integration are to be chosen to make the integral ( p.2| ) well-defined. From 
equation (|3.3|) we have 



TiT = t-t. 



t - 



T T 

T 



t. 



Shifting the variable t by T*T/t, integral ( p.6| ) reads 

It ^ 



dt dt dr* dr 



-{t+T*T/t-t) 



(3.8) 



(3.9) 



The integration over the antisymmetric variables is straightforward, and we are left 
with the t, t integrations. 



1 

2^ 



dt dt 



- t' 



-(t-t) 



(3.10) 



Let us assume for a moment that L is non-integer. Then we may define an 
integration path starting at e^'^cxD, encircling the origin counterclockwise, and returning 
to e^'^oo. The phase angle has to be chosen in a way that the integral is well-defined. 
We end thus with the well-known integral representation for the reciprocal Gamma 
function, both for the t and the t integration, with the result 



27r 



AttL 



r(L + i)r(-L) 

i (e^-^ - l) . 



2sin7r(L + 1) e 



\ttL 



(3.11) 



Equation ( |3.7| ) thus is well-defined for non-integer L if the paths of integration are 
chosen as described above. For integer L the expression on the right hand side is not 
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defined, but it is easily seen that the limit (non-integer L) (integer L) exists and 
gives 

|A|-^ = y J dt ^dtdr*dr iTl^e-T--^^^), (3.12) 



where the i integration is performed counterclockwise on a circle about the origin. 
Equation (|3.12|) holds for all natural numbers L. 
Applied to equation ( |2.24| ), we have 



Z=^ f d[x] d[t] e-t{E-Ho)x|2.|Lg-T4r(i+iAx)] ^3_^3) 
Now the X integrations can be performed, using definition ( p. 23 ), 



1 

IE - EH - XT 



iv)x 



whence follows 



This may alternatively be written as 



Z=l /d[t]e-^[^(^)l, (3.16) 



F{T) = T - L\nT + J2^n (e - E^l - XT) . (3.17) 



where 



Equations ( ^.15| ) to ( |3.17D constitute our next intermediate result. They allow to 
calculate the averaged density of states of a billiard with randomly distributed scatterers 
in terms of the eigenenergies of the billiard without scatterers. All integrals can be solved 
exactly by means of the residuum method. In the remaining step the limit N ^ oo has 
to be performed. The occurring infinite products diverge as a consequence of the delta- 
like singularities in the potential. But the divergencies can be handled in a standard 
way by a renormalization of the coupling constant (Albeverio and Seba 1991). See also 
Bogomolny et al 2001, where the situation of a single scatterer in a rectangular billiard 
is studied. 

From microwave experiments it is known, but for systems with time-reversal 
symmetry only, that in billiards with randomly distributed scatterers the wave functions 
are localized at low energies, but become delocahzed at sufficiently high energies 
(Kudrolli et al 1995, Stockmann et al 2001). A calculation of the averaged density 



of states as a function of energy from equation (3.15), and of two-point correlation 
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function, inverse participation ratio etc. from its generalization should thus exhibit 
clear fingerprints of the localization-delocalization transition. 

This will be the program for future works. For the moment let us proceed along a 
more convenient route by taking fixed and finite, and by performing a second average 
over the shape of the billiard. 



4. The average over the billiard shape 

According to a conjecture of Bohigas, Giannoni, Schmit (1984) the spectrum of a billiard 
with broken time-reversal symmetry should obey the same statistical features as the 
spectrum of a random matrix taken from the Gaussian Unitary Ensemble (GUE). Taken 
this for granted we may replace Hq in equation (|3.13|) by a GUE matrix and perform a 



Gaussian average over the matrix elements to obtain the average over the billiard shape. 
(Up to now Hq had been assumed to be diagonal, but because of the basis invariance 
of the expression we may take any other basis as well; it is much easier to perform the 
average over the matrix elements than over the eigenvalues.) 

The Gaussian average over the matrix elements is trivial and yields 



g-ixtHox^ = /e "'^ ) (4.1) 

/ -i((Ho)2) x;«^,3+?sc;3)(^;:r„+c;cc<)\ 

= (e ) 
= e 2^ ^ ^ . 

Following common practise again we have shifted the average energy to zero, and 
have applied the normalization <^(ifo)^) = N/tt'^ yielding a mean density of states of 
one ed. E = (Verbaarschot et al 1985). After a subsequent Hubbard-Stratonovich 
transformation equation (|4.1|) reads 



where 



y= M^ M . (4.3) 

\T1 y 



To make expression (|4.2| ) well-defined, the y integration has to be performed from 
—00 to 00, and the y integration from — icxo to icxD. Inserting expression ( [4.2|) into 
equation ( |3.13| ) we get as the result of the shape averaging 



1 

h 
1 

h 



d[x]d[t]d[y] iTl^'e-T'-^ 



d[t]d[y] |T|S- 



. TrT -—Try 2 



e - W Ti-^' gixt(E-ATljv)x g-iTr(xy) ^^^^^ 

d[a;]e'"'(^-^^-^)i^". 
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where we have used Tt{XY) = x^Fl^vx. The X integrations are straightforward and 
yield 

(Z) = j-J md[y] mS-^-^e-^^-^^ ^^_^^_^^^ . (4.5) 
Again we apply expression ( |3.7|) to rewrite the determinant, 

^ [d[s] |5|^eT^^(^-^^-^), (4.6) 



E-XT-Y 



^ -TV 



where 



S= !M. (4.7) 



^ a s ^ 

In addition, we replace T by NT, and Y by NY and obtain 

(Z) = -V / m d[y] d[s] \T\'^^-^TrT^-N4TrY- |^|iV ^Tr SiE-XNT-NY) ^4 

= ^ I d[s] \Sf e^'^'^^ y I d[t] |T|^ e-^T^^(i+^^) | d^ 



-N[ ^Try2+Tr(5y) 



The T and the Y integrations can now be performed with the result 



{Z) = ^ f d[g]e^(^^) l*^' ^ e^ '''^ . (4.9) 



|1 + A5|^ 

Equation ( |4.9| ) is the main result of this paper. It is surprisingly simple and allows 
an easy calculation of the density of states of the billiard with randomly scatterers, 
doubly averaged over disorder and shape of the billiard. We only have to perform the 
remaining integrations over 4 commuting and anti-commuting variables. 

5. The density of states 

The calculation of the integral is easiest, if we transform the matrix 5* into a diagonal 
matrix via 

[cr -s J [ -7 VTT^ )[ ■ )[ 7 VI + tP )' ^ ' 

After performing the matrix multiplications we have 

s = sb + {sb - sf) , a* = (3 {sb - sp) , 

o- = -7 {sb - Sf) , s = sf + Pj {sb - sf 

whence follows for the volume element 



(5.2) 



d[s] = -^^^^^^^ . (5.3) 
(sb - Sf) 
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Remembering that the sb integration is from to e^'^oo with a suitably chosen 
phase angle (j), and that the sp integration is along a circle about the origin (see equation 
( P.12| )), we obtain from equation ( |4.9| ) 



X 





f 


i J 






SB 


N 


sf 





dsBdsp d/3 d7 ^j;_^^^^^ff^(^^_,^;)|_^, 



[SB - Sf 



(5.4) 



1 + \Sf 



1 + Xs 



B 



The value for the phase angle can be inferred from equation ( [4.61 ): since t and y are 
real, and has an infinitesimally small positive imaginary part, the integration has to 
be performed from to ioo. 

The integral over the antisymmetric variables is easily done and yields 



1 
27ri 

X 



E+- E'^ 
ip — e 



dSR idsF— -t ^E+SB-E'^sp 



SB 



Sf 



N 



1 + Xsp 



1 + Xsb 



Sb - Sf 

L 

e2^\ B ^F) 



(5.5) 



It follows for the mean density of states (see equation (^ 
{p{E)) = 



1 , dZi 
— Im- 



TT dE' 



(5.6) 



E'=E 



lOO 7-1/ \ 

— Im / dsB f dsp 

vr zvri J J Sb — Sp 





SB 



Sf 



N 



1 + Xsp 



1 + As 



B 



Differentiating with respect to E, we have after some straightforward 
transformations 

(5.7) 



TT 



{p'{E)) = ^lNL{e,a)I(^N-i)L{e,a) 



where 



TT 



/2N 



a 



N/2 



nX 



and 



iMhie, a) 



1 

vri 



dx e^'^'- 



(2x)^ 



In lie, a) = ^ i dye' 

TTl J 



2ey 



[x + a) 
{y + af 



(22/) 



N+l 



(5.8) 

(5.9) 
(5.10) 
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(5.11) 



From the definitions we immediately obtain the recursion relations 

I'nl = I{N+1)L , (-^TVLe^"") = /Ar(L_i)e^^° 

I'nl = —I{N-1)L , (-^AfLe"^^") = — /Ar(L+i)e~^^" 

where the prime denotes differentiation with respect to e. 
For L = we get in particular 



Imie, a) = ^-L- e-'"HN{e) , (5.12) 

^/7^ 



Im{e,a)=^-^H^{e) , (5.13) 

where integral representations of the Hermite polynomials have been used (see e. g. 
Magnus et al 1966). Using the recursion relations we now calculate {p{E)) from equation 
( p. 71) by repeated partial integration with the result 

N-l 

{p{E)) = ^ E hde, a)hL{e, a) . (5.14) 

We have thus obtained a closed expression for the averaged density of states for 
arbitrary values of N and L. It is an easy matter to show, again using the recursion 

oo 

relations ( |5.11| ), that / {p{E)) dE = N, as it should be. 



For L = equation (|5.14|) reduces to 



{p{E)) = ^ E [Me)]' , (5.15) 



where 



is an harmonic oscillator eigenfunction. This is identical with the well-known exact 
expression for the density of states of the Gaussian unitary ensemble, which in the limit 
of large N reduces to Wigner's semicircle law (Mehta 1991). 

6. The strong coupling limit 



Using the recursion relations (|5.11|) Inl and Inl can be calculated from Ijsfo and I no 
by repeated integration or differentiation, respectively. Since all integrations can be 
performed analytically, we have got an exact representation for the density of states for 
arbitrary L. Though this may be helpful for small values of L, it is not very useful for 
practical purposes, since one is usually interested in the limit N, L ^ oo while the ratio 
I = L/N remains finite. 

In such a situation it suggests itself to solve the integrals (|5.9|) and (|5.1CI| ) with 



help of saddle point techniques. This leads to a cubic saddle-point equation which 
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still can be solved exactly using Cardano's formula. The resulting equations are not 
very elucidating, however. Therefore we proceed in another direction and restrict the 
following discussion to the strong coupling limit A ^ 1, or a -C 1. In the discussion we 
have to discriminate between the two situations N > L and N < L. 
(i) N> L 

For this case we may replace {x + a)^ and {y + a)^ in the integrands by and y^, 
respectively, to obtain 

/^^(e,a) = _ / dxe2"'(2x)'^"^e"% (6.1) 

TTl J 

— ioo 

lNL{e,a) = — Uye-'^y{2yr^''-'^^e-y\ (6.2) 

TTl J 

In the strong coupling limit the averaged density of states for a billiard system with 
N levels taken into account and L randomly distributed scatterers is thus the same as 
for a system with N — L levels, and no scatterer at all. We are again in the random 
matrix regime. 

Remember that already in the beginning we observed that the spectra of billiards 
with randomly distributed scatterers show random matrix behaviour, but at that point 
we considered the weak coupling limit A <^ 1 ( see the discussion following equation 

(ill))- 



We thus can note that for > L both in the weak and the strong coupling limit 
the averaged density of states shows random matrix behaviour, 
(ii) N <L 

Now we cannot replace any longer (x + a)^ and [y + a)^ in the integrands in 
equations ( |5.9| ) and ( |5.10| ) by x^ and y^, since in this limit the integral for Inl diverges, 
and that for Ij^^ gives zero. In the limit a <^ 1, on the other hand, the main 
contributions to the integrals come from regions a; -C 1 and y <^1, where the Gaussian 

2 2 

cut-offs are not yet relevant. We may therefore replace e^ and e~^ by one, and solve 
the integrals by means of the residuum method with the result 



2-'^ (^j (e-.~)U„ (6.3) 

= 2^+ia^+i-^e(^)-^e-^z^-^-i4^-^-^)(^) , 
[L - 1)! 



L-l 



N 



^ N\ I dy, 



{y + a) 



j/=0 



(6.4) 



-{L-N). 
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where 

z = 2ea = E/\. (6.5) 

Q{z) is the Heaviside step function, and LI^\z) is a generahzed Laguerre 
polynomial. (There are two conventions for the Laguerre polynomials found in literature, 
differing in the normalization. In this paper the definition of Magnus et al 1966 is 
adopted, where L^")(0) = ("^°)-) It follows from equation ( |5.14| ) for the density of 
states 

1 N-l 

(piE)) = T77^0(^) E k\ z'-'-'Lt'-'\z) Lf-'\z) . (6.6) 



Equation ( |6.6|) simplifies considerably in the limit L ^ oo, — oo, with L/N 
remaining finite. Inserting expressions ( |6.3| ) and ( |6.4| ) for lNL{z,a) and lNL{z,a), 
respectively, into equation ( |577| ) we obtain 



(p'iE)) = ^ e(.) e-z--^-Lr^-^)(z) L[i,-Jl,^'\z) . (6.7) 



In terms of the function 

z a + 1 



i^\z) = e-h^L^:\z), (6.8) 



equation (|6.71 ) may be written as 



ip'iE)) = ^ (^e(.).-,r"-^)(z),^_-r^)(z). (6.9) 

The 7/^°-' obey the differential equation 

„ f2n + a + l 1 l-a'^\ 



Equation (|6.10| ) is easily identified as the radial Schrodinger equation of the 



hydrogen atom, where y'^\z)/z is the radial part of the wave function. This suggests 
an approximation of y^'' by means of the WKB method. In the present context it is 
sufficient to consider the solution in the classically allowed region. For this regime the 
WKB approximation yields (see e. g. section 9.3 of Morse and Feshbach 1953) 

y^:\z) = ^cos \ /gdz-;i , (6.11) 

where 




2n + a + l 1 a2 
and zi are the classical turning points given by 



Zq/i = 2n + a + l± J{2n + l){2n + 2a + 1) . (6.13) 
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The replacement of 1 — by —a^ in going from equation ( |6.10| ) to equation ( |6.12| ) 
corrects for the singularity of the potential at z = 0, see the discussion in Morse und 
Feshbach 1953. (The same technique can be applied to derive the semi-circle law in a 
simple way from the exact expression ( |5.15| ), see chapter 3.2.3 of Stockmann 1999; the 
procedure is more or less an elaboration of an idea developed in appendix A. 9 of the 
book of Mehta (1991).) 

Inserting approximation (|6.11|) into equation ( p. 91) we end up with 

1 1 



{p{E)) = —- ^ALN -{z-L- Ny, z = E/\, (6.14) 
Att Iz * 

where the classical turning points are given by 



zq/i =L + N± 2VLN . (6.15) 

Details of the derivation can be found in the appendix. The density of states thus 
changes dramatically if L surpasses N. For L < N Wigner's semicircle law is found, and 
the eigenenergies are distributed between —2tt/N and 2tt/N . For L > N on the other 
hand only positive eigenvalues are found, if A is positive, in an energy window limited 
by A^^o and Xzi. 



7. Discussion 



We have obtained a surprisingly simple expression for the averaged density of states of 
a billiard with randomly distributed scatterers. The central ingredient was the idea to 
substitute the average over the scatterer positions in equation ( p.l6| ) by an weighted 
average, with the wave function amplitude probability density as the weight function. 
It was argued that both averages are equivalent. In view of the central importance of 
this procedure it seems appropriate to discuss the limitations of the approach. 

(i) First, the impurities are considered as uncorrelated. In particular, it is not 
excluded that two impurities occupy the same site. 

(ii) Second, wave functions belonging to different eigenvalues are considered as 
uncorrelated. This may pose a problem, since it is known from semi-classical quantum 
mechanics that there are correlations on energy scales of the order of h/T, where T is 
the length of the shortest periodic orbit (see Gutzwiller 1990 for a review). On the other 
hand, these correlations vanish in the semi-classical limit on energy scales of the mean 
level spacing. It therefore seems legitimate to neglect correlations between different 
wave functions. 

No problem, on the other hand, arises from the fact that there are spatial 
correlations for individual wave functions, as is well-known from the works of Berry 
(1977) and Fal'ko and Efetov (199). Since only the weight of the wave function 
amplitudes enters equation ( ^.18| ), spatial correlations are completely irrelevant. 

The approximation performed in section ^ by substituting the billiard by a random 
matrix of finite rank is of another type. It has been applied to obtain a simple tractable 
model, but by this second step we have reduced our system to a mere caricature of a real 
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billiard system. In particular the information on the dimension of the billiard, which 
obviously is an important quantity for questions of localization and delocalization, is 
lost. (The information on the dimension is still present in equation ( p.l5 ), namely in the 
spectrum of the empty billiard which depends on the dimension via the mean density 
of states.) 

This is why in the moment a comparison with literature results is not possible. In 
particular the work of Luttinger and Tao (1983) has to be mentioned in this respect, 
who calculated the density of states for the billiard with randomly distributed scatterers 
in the low energy limit. For a more detailed consideration of their results, we would 
have to go back to equations ( p.l5[ ) to ( |3.17D , perform the limit cxd, and calculate 
the density of states for the true billiard system, and not a random matrix substitute 
only. 

But the present results suggest that already in our toy model there is a localization- 
delocalization transition at L = A^. For L < N we are in the regime of delocalized wave 
functions obeying random-matrix behaviour. For L > N, on the other hand the wave 
functions become localized, giving rise to a completely changed density of states. In 
this respect the rank of the matrix seems to take the role of the energy in the real 
billiard system. 

For the moment, however, this conclusion must be considered as premature. 
Knowledge of the density of states is not sufficient to discriminate between localized 
and delocalized wave functions. For this we need additional information on the two- 
point correlation function, the inverse participation ratio and related quantities. The 
corresponding studies are under progress and will be published separately (Guhr and 
Stockmann 2002). 

If one compares the present approach with the non-linear a model, a dramatic 
simplification is found. In the non-linear a model one ends up with a supersymmetric 
integral over supersymmetric field variables which can be solved only within the zero- 
mode approximation. In our approach the very simple integral ( [4.9|) is obtained instead, 
containing only one set of supersymmetric variables, which for the density of states even 
can be solved exactly. The same is true for all ra-point correlation functions as will be 
shown in Guhr and Stockmann (2002). 

Even better, the assumption of a random distribution of point-like scatterers applied 
in this work is a much more realistic description of the situation found in mesoscopic 
systems than the assumption of a delta-correlated disorder potential assumed in the 
non-linear a model. 

It might be considered as a draw-back that the present derivation is based on two 
unproven conjectures, namely that (i) the wave function amplitudes in a chaotic billiard 
are Gaussian distributed, and that (ii) the eigenvalues in a chaotic billiard obey random 
matrix behaviour. On the other hand, there is such an overwhelming numerical evidence 
that both conjectures are true that one could equally well argue that both assumptions 
are even better founded than the assumption of Gaussian distributed matrix elements 
applied in random matrix theory. 
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From the point of view of an experimentalist it would be highly desirable, if all 
quantities of interest are available for systems with time-reversal symmetry as well. 
Though the calculations for this case are notoriously difficult, it should be worthwhile 
to undertake the effort. Experiments with microwaves on localization-delocalization 
transitions, pulse propagation etc. in disordered systems do already exist, as was 
mentioned in the introduction, and wait for their proper theoretical explanation. 
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Appendix A. The density of states in the large L limit 

To derive equation ( |6.14| ) for the density of states we start with equation (|6.7|), 



e z 



L-N-l j{L^N-l) 



z]L 



(L-N+l), 

N-1 y 



(A.l) 



where we have introduced z = E/X as a new variable, and where p{z) dz = p{E) dE. z 
is assumed to be positive in the following. Using elementary relations for the Laguerre 
polynomials, equation ([A.l|) may be transformed as follows: 



N\ 
N\ 



L-N jj^L-N) 



N 



-(L-N) 



62 z 2 



{L-N) 



L-N+l 



62 Z 



Vn 



{L-N) 



(A.2) 



where is given by equation ( |6.8|) . It remains to determine the normalization 

constant yo. From the orthogonality relation for the Laguerre polynomials we have 

2 dz 



6 ^X" 



dz 



[n + a)\ 



From the WKB approximation (|6.11|) , on the other hand, we obtain 



2 dz _yl /-^i dz 
z 2 Jzo zq{z) 



(A.3) 



(A.4) 



whence follows 



I/O 



1 + a)! 



(A.5) 
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The sign has to be chosen positive to be in accordance with the usual definition 
of the Laguerre polynomials. Inserting now the WKB approximation for y^^\z) into 
equation (|A.2|) we have 



p'{z) = —e 
n 



z_ L-N-l 1 

Q2 z 2 cosw 



z L-N+1 1 

62 Z 2 COSW 



where 



and 



q{z) 



lL + N + 1 _1 _{L-Ny 



2z 



4^2 



w 



/ q{z) dz--. 



It follows 



. 1 \fl , L-N-l lq'\ 1 " 

p (Z) = H -— cosw — vosmw 

^ TTZ [\2 2z 2q J y/q 

'1 L-N + 1 lq'\ 1 

— cos w — a/o sm w 

2 2z 2q) ^ ^ . 



X 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



The q' terms may be discarded since they are by an order of 1/L smaller as the 
other ones. For the same reason we may replace L-N-l and L — N + lhyL — N. 
Averaging over the rapidly oscillating terms, we have 



27izq 



1 L-N\ n L-N 

2 ^ 2z ) \2 2l~ 



(A.IO) 



Inserting expression ( |A.7] ) for q, we obtain 



2TTzq 

irn ) 



L + N {L-Nf 



2'Kq 



2z 

= -q ■ 

IT 



2z^ 



where again terms of the order of 1/L have been neglected. It follows 
1 



TC 



which is equivalent with equation (|6.14| ). q. e. d. 



(A.ii; 



(A.12) 
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